A&A manuscript no. 

(will be inserted by hand later] 



Your thesaurus codes are: 

05(10.07.3 uCen, NGC6397, NGC6752, Lillerl, 13.25.5) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Multiple and variable X-ray sources in the globular clusters 
^ ,u;Cen, NGC 6397, NGC6752, and Liller 1 

F. Verbunt^ and H.M. Johnston^ 



o 
o 



Astronomical Institute, P.O.Box 80000, NL-3508 TA Utrecht, The Netherlands 



O ^ ^ Anglo Australian Observatory, P.O. Box 296, Epping, NSW 2121, Australia 



> 

o ' 

o . 
o . 

o ■ 

Oh; 

6 ■ 



X 



February 5, 2008, accepted 



Abstract. We detect stars from the Hipparcos and Ty- 
cho Catalogues in the field of view of observations with 
the ROSAT HRI of three globular clusters. We use the 
positions of these stars to reduce the systematic error in 
the positions of X-ray sources in the clusters to ~ 2" for 
wCen and NGC 6752, and 1" for NGC 6397. We detect 
three X-ray sources in the core of u> Can, and four in the 
core of NGC 6752; the data for the center of NGC 6397 
may be fitted with five or six sources. Outside the cores, 
but within the half-mass radius of the clusters, we de- 
tect two sources in uj Cen, one in NGC 6397 and two in 
NGC 6752; these may or may not be cluster members. A 
ROSAT HRI observation of Liller 1 does not detect a low- 
luminosity source, at a limit below a detection with ASCA. 
We discuss the nature of the low-luminosity X-ray sources 
in globular clusters in the light of these new results. 

Key words: globular clusters: individual oj Cen, 
NGC 6397, NGC 6752, Lillerl - X-rays: stars 



1. Introduction 

Globular clusters contain many X-ray sources at lower 
luminosities, ^ 10^'*ergs~^. These sources were first 
discovered with the Einstein satellite (Hertz & Grindlay 
1983), and many more were found with ROSAT (for a 
compilation, see Johnston & Verbunt 1996). The nature 
of these low-luminosity sources is the subject of debate, 
because various types of objects can emit X-rays at such 
luminosities, such as soft X-ray transients in quiescence, 
cataclysmic variables, RS CVn binaries, and recycled neu- 
tron stars (see e.g. Fig. 8 in Verbunt et al. 1997). The 
most compelling identification of a dim X-ray source with 
an object observed at other wavelengths is the recycled 
radio pulsar in M28: the X-ray fiux varies on the pulse 
period (Banner et al. 1994). Plausible identifications with 
cataclysmic variables have been suggested for dim X-ray 
sources in NGC 6397, NGC 6752, NGC 5904 and 47Tuc 



(Cool et al. 1995b, Grindlay 1993, Hakala et al. 1997, Ver- 
bunt & Hasinger 1998). These identifications are based 
on the proximity of the X-ray position to that of a cat- 
aclysmic variable, and thus their probability depends on 
the accuracy of the X-ray position. 

The accuracy of the ROSAT position of a detected 
X-ray source is determined by two factors: the statisti- 
cal accuracy of the position of the source on the detector, 
and the accuracy with which the position of the detector 
as a whole is projected on the sky. For a sufficient number 
of photons the statistical error is less than an arcsecond, 
but the projection in general has a typical error of 5". 
Secure identification of a source in the detector field re- 
duces the error in the projection to the statistical error 
of the identified source, provided that the optical (or ra- 
dio) position has better accuracy. Only one identification 
is necessary, because the roll angle of the detector (i.e. 
the North-South direction) is accurately known; nonethe- 
less, identification of more than one source is preferable to 
allow checks on internal consistency. In a globular cluster 
the surface density of possible counterparts is so high that 
chance coincidence usually cannot be excluded; a secure 
identification can usually be made only for X-ray sources 
detected well outside the cluster. This method has been 
used by Verbunt & Hasinger (1998) to improve the posi- 
tional accuracy of the sources in the core of 47 Tuc from 5" 
to 2", whereby the area in which the source is expected 
to lie is reduced sufficiently to exclude several proposed 
identifications, and increase the probability of others, in- 
cluding two possible cataclysmic variables. 

In this paper we investigate three clusters known to 
contain multiple dim X-ray sources in their core, which 
have been observed in long exposures with the ROSAT 
HRI, and one cluster known to harbour a transient. We 
analyse hitherto unpublished observations and detect both 
previously published and new X-ray sources. All source 
positions are checked in the SIMBAD data base versus po- 
sitions of other objects, and we find objects in the Hippar- 
cos or Tycho Catalogues (ESA 1997, Ferryman et al. 1997, 
H0g et al. 1997) with each cluster, i.e. counterparts with 
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Table 1. Log of the ROSAT HRI observations of globular 
clusters analysed in this paper. For each cluster observa- 
tion, the observation date(s) and exposure time are given. 
We further give the shift in a, i5 applied to bring the X-ray 
coordinate frame of the longest observation to the optical 
coordinate frame J2000. 
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12958 






1995 


Mar 
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NGC 6752 
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1995 
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1996 
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2450182.957-186.041 
Liller 1 


23278 






1996 


Sep 


2450332.541-337.810 


16412 







describe the resuhs for uj Cen, NGC 6397, and NGC 6752, 
respectively. In Sect. 6 we discuss an observation of Liller 1. 
A discussion of our results is given in Sect. 7. 

2. Observations and data reduction 

The X-ray observations were obtained with the ROSAT X- 
ray telescope (Triimper et al. 1991) in combination with 
the high-resolution imager (HRI, David et al. 1995). The 
list of the observations is given in Table [l[ The stan- 
dard data reduction was done with the Extended Scientific 
Analysis System (Zimmermann et al. 1996), as follows. 
To take into account the re-calibration of the pixel size 
(Hasinger et al. 1998), we multiply the x,y pixel coordi- 
nates of each photon with respect to the HRI center with 
0.9972. A search for sources is made by comparing counts 
in a box with the counts in a ring surrounding it, and by 
moving this detection box across the image. The sources 
thus detected are excised from the image and a back- 
ground map is made for the remaining photons. A search 
for sources is then made by comparing the number of pho- 
tons in a moving box with respect to the number expected 
on the basis of the background map. Finally, at each po- 
sition in which a source was found, a maximum-likelihood 
technique is used to compare the observed photon distri- 
bution with the point spread function of the HRI (Crud- 
dace et al. 1988). This produces a maximum-likelihood 
value ML such that the probability that the source is due 



sure that all such sources are found, we enter in the max- 
imum likelihood technique all sources that have ML>10 
according to the sliding box searches.) 

Upper limits at the position of known sources were de- 
termined by counting the number n of actually detected 
photons at the position of the source (and in an area sur- 
rounding it corresponding to the uncertainty in the posi- 
tion); we then assign as upper limit the lowest expected 
number m for which the probability of measuring a num- 
ber n or smaller is less than 5% according to the Poisson 
distribution. 

The maximum-likelihood technique also provides an 
indication whether the source is extended. If such indica- 
tion is present, we apply further analysis to test whether 
the source is a multiple point source. 

The further analysis is also based on maximum like- 
lihood techniques, but the analysis is limited to a small 
area of the detector, near its center. This allows the sim- 
plifications that the background in the analysed area is 
a constant (as opposed to a polynomial function of the 
X, y pixel coordinates), and that the point spread function 
is that for the center of the image (David et al. 1995). 
Suppose that a model to be tested predicts photons 
at detector pixel i. The probability that Ui photons are 
observed is then given by the Poisson probability: 

P. ^ ^^^^ (1) 

The probability that the model describes the observations 
is given by the product of the probabilities for all i in the 
region considered: L' — WPi. For computational ease we 
maximize the logarithm of this quantity: 

In L' = In Pi ~ rii In m; — — In n;! (2) 

i i i i 

The last term in this equation doesn't depend on the as- 
sumed model, and - in terms of selecting the best model 
- may be considered as a constant. Thus maximizing L' 
is equivalent to minimizing L, where 

ln_L = — 2(^~^ fij hrnij — mj) (3) 

i i 

If one compares two models A and B, with number of fitted 
parameters ua and n b and with likelihoods of In La and 
In Lb, respectively, the difference AL = Ini^ — In is is 
a distribution with ua — ns degrees of freedom, for 
a sufficient number of photons (Cash 1979, Mattox et al. 
1996). 

Our analysis of possibly multiple sources thus proceeds 
as follows. First we compute AL for a model with constant 
background and for the best model with background plus 
one source, and compare it with the x^-distribution with 3 
degrees of freedom. If AL > 15, the presence of one source 
has a significance more than three sigma. Next we compare 
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Fig. 1. X-ray contours in the central area of uj Cen as ob- 
served with the ROSAT HRI in 1995 July. The image was 
smoothed with a 2-d cr~5" Gaussian. The detected sources 
are indicated with their numbers in Table |^. The inner cir- 
cle gives the core radius of the cluster, the outer circle the 
half-mass radius. The lower and left axes give pixel num- 
bers for the ROSAT HRI detector, the upper and right 
axes right ascension and declination with respect to the 
cluster center. The conversion between pixel and celestial 
coordinates is accurate to within 2". 

the best models with three and two sources to prove the 
significance of a third source, etc. until no more significant 
sources are found. 

The addition of one source adds three fitted parame- 
ters, one for its number of counts and two for its position. 
In the case of NGC 6397 optical counterparts have been 
suggested for three X-ray sources. For these we also make 
a fit in which the distances in right ascension and declina- 
tion between these three sources is fixed to the optically 
determined values. The three sources in that case only add 
five fitted parameters, two for the position of one of them, 
and three for the fluxes. 

To determine the error in a parameter, we start from 
the best fit value a^. We then fix the parameter at -|- d 
and make a new fit, allowing all other parameters to vary. 
The value of d for which In L increases by 1 is quoted as 
the 1-sigma error. 

3. LJ Cen 

Lu Cen is a very massive globular cluster, with a relatively 



Table 2. X-ray sources detected in the globular cluster 
LoCen {Av = 0.47, d = 4.9 kpc, Djorgovski 1993) with 
the ROSAT HRI. Sources with a number less than 17 cor- 
respond to sources detected previously with the ROSAT 
PSPC (Johnston et al. 1994); those with a higher num- 
ber are new sources. For each source we give the position, 
the statistical error in the position (in "), the distance 
to the cluster center in units of the core radius rc, the 
countrate with error, and where applicable the identifica- 
tion with sources detected with Einstein (Hertz & Grind- 
lay 1983). The sources are ordered on declination. The 
positions given are those after correction for bore sight 
(see text). The positions of the center of the cluster (GC, 
Djorgovski & Meylan 1993), its core radius and half-mass 
radius (Trager et al. 1993), and positions of some optical 
objects discussed in the text are also listed; epochs are 
1990.5 for positions by Cool et al., 1979 for USN0-A2, 
and 1995.5 for HD 116789. 



X-ray sources 



X 


Q 


(2000) 


S (2000) 


A 


d/rc 


cts/ksec 


4 


13 


27 27.72 


-47 


19 


8.0 


0.4 


4.6 


7.8±0.4 D 


3 


13 


25 52.05 


-47 


19 


8.6 


0.5 


5.1 


6.3±0.3 A 


18 


13 


26 46.33 


-47 


19 


45.9 


1.0 


3.4 


0.5±0.1 


19 


13 


27 21.21 


-47 


23 


22.0 


1.8 


3.1 


0.4±0.1 


6 


13 


27 29.27 


-47 


25 


54.5 


1.3 


3.0 


0.9±0.2 


20 


13 


26 48.70 


-47 


27 


46.5 


1.0 


0.4 


0.6±0.1 


9a 


13 


26 53.39 


-47 


29 


1.5 


0.9 


0.5 


0.9±0.1 C 


7 


13 


26 19.75 


-47 


29 


11.9 


0.8 


1.7 


1.4±0.2 B 


9b 


13 


26 52.04 


-47 


29 


37.2 


0.9 


0.6 


1.1±0.2 C 


21 


13 


26 48.99 


-47 


31 


28.1 


2.1 


1.1 


0.5±0.1 


10 


13 


26 25.20 


-47 


32 


28.9 


1.4 


2.0 


0.7±0.1 


13 


13 


26 11.12 


-47 


37 


11.1 


3.7 


4.0 


1.3±0.3 






optical objects 












GC 


13 


26 45.9 


-47 


28 


37 




155", 


Th = 288" 


4 


13 


27 27.37 


-47 


19 


6.2 


Cool et al. (1995a) 


4 


13 


27 27.68 


-47 


19 


6.3 


USNO-A2 0375-18249604 


3 


13 


25 51.79 


-47 


19 


7.0 


Cool et al. (1995a) 


3 


13 


25 52.23 


-47 


19 


7.5 


USNO-A2 0375-18177834 


18 


13 


26 46.33 


-47 


19 


45.9 


HD 116789 




E 


13 


29 18.54 


-47 


22 


50.5 


Marf 


ron & Bolte (1987) 



outside the cluster core, and it appears that only source C 
is clearly related to the globular cluster (e.g. Verbunt et 
al. 1995). Source E has tentatively been identified with a 
foreground K star (Margon & Bolte 1987). The sources A 
and D have been identified with foreground M stars, on 
the basis of better positions for the X-ray sources obtained 
with Einstein and ROSAT HRI observations (Cool et al. 
1995a). A ROSAT PSPC pointing indicates that source C, 
near the cluster center, is composed of two sources of com- 
parable luminosity (Johnston et al. 1994). 

We have analysed all observations listed in Table ^ 
The 1992 and 1993 data have been reported on before by 
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1995. We use this observation as the basis for comparison 
with the other observations. 

3.1. Source list and membership probability 

In the July 1995 observation we detect twelve sources, 
listed in Table ^ Eight of these sources have been detected 
before with Einstein or with the ROSAT PSPC observa- 
tion reported by Johnston et al. (1994); four sources are 
new. 

The sources X 3 and X 4 arc the Einstein sources A and 
D, respectively, identified with foreground M dwarfs by 
Cool et al. (1995a, see Table|^). Both stars can be found in 
the USN0-A2 catalogue (Monet et al. 1998, see Table |). 
The new source X18 can be identified with HD 116789; 
this is star TYC 8252 46271 in the Tycho Reference Cat- 
alogue, and thus its position and proper motion are well 
known (H0g et al. 1998). For this reason we use this star 
to determine the bore sight correction, i.e. the offset be- 
tween the X-ray coordinates and the optical coordinates. 
The result is listed in Tabled. 

This bore sight correction has been applied to the X- 
ray positions, and the resulting positions are given in Ta- 
ble^ With a statistical accuracy for the X-ray position 
of X 18 of ~ 1", and taking into account small additional 
systematic errors (see Hasinger et al. 1998), we estimate 
the systematic error in the positions given in Table || to 
be less than 2"; for each individual source its statistical 
uncertainty should be added in quadrature to this sys- 
tematic error. The positional accuracy can be improved 
if accurate astrometry of X3/A and X4/D is obtained, 
which will allow computation of their positions at epoch 
1995.5. 

In the ROSAT Deep Survey (Hasinger et al. 1998), an 
area with 12.'5 radius contains 25 sources brighter than our 
approximate detection limit of 0.7 x 10^^'* ergcm^^s^^; 
we thus expect '^1 serendipitous source within the core 
radius, i.e. the faintest source within the core radius, 
X20, may well be a serendipitous background source. The 
brighter sources within the core radius probably are as- 
sociated with w Cen. Outside the core radius a source is 
more likely to be a fore- or background source than a clus- 
ter member. 

3.2. Sources in the cluster 

The X-ray image of the inner area of w Cen is shown in 
Fig.|l|. The half- mass radius of the cluster contains five 
sources. Source X9 from Johnston et al. (1994) is clearly 
separated into two sources, which we denote X 9a and X 9b 
for the northern and southern source, respectively. 

We have determined the countrates or upper limits for 
the six central sources in all ROSAT HRI observations of 
Lo Cen listed in TablepI For the very short observations, no 



four of the central sources in Fig.g. We also show the 
PSPC observation, dividing the PSPC counts for PSPC- 
X9 equally between X9a and X9b. For the absorption 
towards w Cen and a black body spectrum of 0.6 keV the 
PSPC countrate is about 2.8 times the HRI countrate; 
for smaller reddening the PSPC-to-HRI count ratio varies 
rapidly, and therefore we do not show PSPC points of the 
foreground M dwarfs, whose absorption is unknown. 

There is marginal evidence for variability in sources 
X9b and XIO; and no evidence for variability of X7 and 
X9a. X20 is detected only in July 1995 and in 1996, X21 
only in July 1995, and all upper limits in the other ob- 
servations are compatible with the faint fluxes of these 
sources. 

The July 1995 observation was obtained in two parts, 
separated by about 10 days. We have analysed the two 
parts separately, and find marginal variation of X7 (1.7 ± 
0.3 and 1.1 ±0.2 ctsksec"^ in the first and second half, re- 
spectively) and of X9b (0.8 ±0.2 and 1.4 ±0.2 ctsksec"!, 
respectively). In addition, virtually all counts of X21 are 
from the first part of the observation. The number of 
counts of X21 is too small for further subdivision. 

To convert the observed countrates into X-ray lumi- 
nosities, we assume a column and distance to oj Cen as 
given m Table |. For an 0.6 keV blackbody (see the anal- 
ysis of the PSPC spectrum of X9 in Johnston et al. 
1994) 1 ctsksec"^ in the ROSAT HRI corresponds to 
1.5 X 10'^^ ergs^^ in the 0.5-2.5 keV band. The two sources 
X 9a and X 9b thus have about this luminosity, source X 7 
is 40% brighter, and sources X 20 and X 21 are 40% fainter. 

3.3. Sources not related to the cluster 

The foreground dwarfs (Einstein X-ray sources A and D) 
are highly variable, as has been pointed out before (Koch- 
Miramond & Auriere 1987, Cool et al. 1995a). These 
sources also vary between the first and second half of the 
July 1995 observation. The extremely high flux of A in 
1996 is due to a flare, which lasts almost a day (see Fig. ||). 

HD 116789 is an AOV star. From its magnitude V = 
8.40 and colour B-V = 0.07 we estimate Ay = 0.22 and a 
distance of about 310 pc. For an assumed bremsstrahlung 
spectrum of 1.4 keV the observed countrate corresponds to 
a luminosity in the 0.5-2.5 keV band of ~ 2 x 10^^ ergs^^. 
It is not expected for an AOV star to emit such a flux; 
perhaps this star has a white dwarf companion which 
is responsible for the X-ray emission, as various other 
AOV stars detected with ROSAT; on the other hand, var- 
ious apparently single AOV stars in the Bright Star Cata- 
logue have been detected at similar and higher luminosi- 
ties as HD 116789 (e.g. HD 17864, Hunsch ct al. 1998). 
HD 116789 was detected with EXOSAT by Verbunt et al. 
(1986), who interpreted the detection as due to the ul- 
traviolet leak of the CMA detector. The EXOSAT CMA 
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Fig. 2. X-ray lightcurvcs of four sources in the center of u; Cen and of two foreground M dwarfs (A, D). The HRI data 
are shown as • (detections, with 1-cr errors), | (2-0- upper limits). PSPC data (converted to estimated HRI countrates) 
are shown as o. There is marginal evidence for variability in sources X9b and XIO; the foreground sources A and D 
are highly variable, with source A showing a large flare in 1996. 



~ 10'^"'^ erg s""'^, much higher than the luminosity derived 
for this star from the ROSAT observations; we therefore 
still think that the EXOSAT countrate is due to the ultra- 
violet flux. The ultraviolet leak in the ROSAT HRI is far 
too small (Berghofer et al. 1999) to explain the ROSAT 
detection. 



o 
CO 



o 



o 
o 



o 



3.4- Discussion 

The positions of the sources within the half-mass radius 
of uj Cen as given in Table I are more accurate than previ- 
ously published positions, and may be used to search for 
optical counterparts. We have done this among the vari- 
ables (contact binaries, detached binaries, and suspected 
RS CVn stars) found in w Cen by Kaluzny et al. (1996, 
1997): no counterpart is among these stars. (Only one 
of these variables is in the area shown in Fig.|l|, viz. the 
contact binary 0GLEGC13.) Our non-detection of these 
binaries is not surprising, considering that our detection 
limit is above 10^^ ergs~^: all of the contact binaries hith- 
erto detected in X-rays (McGale et al. 1996), and many 



o 



Jijt 



129.5 130 
time (JD-2450000) 



1 30.5 



Fig. 3. Lightcurve of X3/A in the 1996 observation. • 
indicate the numbers of counts collected in 100 s intervals; 
the horizontal lines indicate when ROSAT was collecting 
data. 



Cool et al. (1995a) argue that X7/B is an extended 
source. This source is detected in the ROSAT PSPC ob- 
servation (Johnston et al. 1994) and in the ROSAT HRI 
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and 1993 observations used by Cool et al. (1995a); in all 
of them X7/B is compatible with being a point source. 

With the identification of sources X3 and X4 with 
foreground stars, we can reinvestigate the suggested iden- 
tification of X5/E with a foreground K star, as suggested 
by Margon & Bolte (1987). We use the optical positions of 
A and D to determine the offset between the X-ray posi- 
tions of the PSPC observation as listed in Johnston et al. 
(1994) and the optical positions. We then apply this offset 
to the position of X 5, and find that the resulting position 
is at 6±5 arcseconds from the optical star. We take the 
position of the optical star (Table|) from USNO A2 0375- 
18334783 (Monet et al. 1998). Identification of X5 with 
the foreground star is therefore a distinct possibility. 



4. NGC6397 

NGC6397 is a nearby cluster, with a collapsed core, in 
or close to which Cool et al. (1993) detected four X-ray 
sources (B, Cl-3) with a ROSAT HRI observation. Pho- 
tometry with the Hubble Space Telescope enabled Cool et 
al. (1995b) to find eight candidate counterparts for these 
sources, on the basis of high ultraviolet flux or of H a emis- 
sion. The Ha emission of three stars has been confirmed 
spectroscopically by Grindlay et al. (1995) who argue that 
these stars are cataclysmic variables, and responsible for 
the X-ray emission close to the core. 

4-.1. Source list and membership 

We analyse first the longest observation, obtained in 1995, 
and use this as a reference for our discussion of the ear- 
lier, shorter observations. The standard analysis provides 
14 sources, listed in Table ^. Identifications with earlier 
X-ray sources or optical objects are indicated; 7 sources 
are new. X6 has been identified by Cool et al. (1993) as 
SAO 244944. This star is identical to HD 160177, and is in 
the Hipparcos Catalogue as HIP 86569. Its position and 
proper motion are thus very accurately known, and we use 
it to determine the bore sight correction. This bore sight 
correction is given in Table|l|, and is applied to the X-ray 
positions; the resulting positions are given in Table |. The 
statistical uncertainty in the X-ray position of X 6 is about 
0.5"; we therefore estimate that systematic error of the X- 
ray positions listed in Table ^ is better than 1"; this error 
should be added in quadrature to the statistical error for 
each individual source position. The quasar identified by 
Cool et al. (1993) with X5 coincides within the error with 
our position for X5. However, the active galaxy identified 
by Cool et al. (1993) with X2 is 10" from our X-ray po- 
sition, mainly in right ascension; and we conclude that it 
is not the X-ray source. The explanation probably lies in 
the new scale for the size of the HRI pixel that we use 



Table 3. X-ray sources detected in the NGC6397 {Ay = 
0.56, d = 2.2 kpc, Djorgovski 1993) with the ROSAT HRI, 
for the standard analysis of the whole field, and separately 
for two multi-source analyses of the central area. Numbers 
up to 10 are sources from Johnston et al. (1994), higher 
numbers are new; cross-identifications with sources dis- 
cussed by Cool et al. (1993) are listed on the right. All 
X-ray positions have been corrected for boresight. The 
positions of the center of the cluster (GC, Djorgovski & 
Meylan 1993), its core radius and half-mass radius (Trager 
et al. 1993) and the positions of some optical objects dis- 
cussed in the text are also listed; epochs are 1992.7 for 
positions by Cool et al., and 1996.3 for HIP 86569. 
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The flux detection limit is about 0.8 x 
10^^^ ergcm~^s~^ outside the blended central re- 
gion, similar to that obtained for uj Cen. Analogous to 
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xpix 

Fig. 4. Positions of ultraviolet or H a-emission objects in 
the central area of NGC 6397 (x , +, numbered with their 
ID in Table 2 of Cool et al. 1995b) superposed on X-ray 
contours of sources X13/B and X4/C as observed with 
the ROSAT HRI in 1995. The candidate counterparts for 
three X-ray sources suggested by Cool et al. are marked 
X . The X-ray image was smoothed with a 2-d cr~2" Gaus- 
sian. The lower and left axes give pixel numbers for the 
ROSAT HRI detector, the upper and right axes right as- 
cension and declination with respect to the cluster center. 
The conversion between pixel and celestial coordinates is 
accurate to within 1". 



expect 1.4 background sources within 3' from the center 
of of NGC 6397; the sources at 0.'5 < r < 3' therefore 
may be background sources. We thus cannot decide 
whether X 12 is a cluster member. Outside the half-mass 
radius, the sources are more likely to be background or 
foreground sources. X5, just outside the half-mass radius, 
is a quasar (Cool et al. 1993). 

4-2. The central sources 

In Fig.H we show the X-ray contours of the center 
of NGC 6397 together with the ultraviolet and/or Ha- 
emission stars discovered by Cool et al. (1995b). The 
first models we investigated as fits to the central area of 
60" X 60", containing sources X 13/B and X4/C, are those 
with successively one, two, three, four and five sources; all 
with free positions. Using the AL criterion for significance 
(see Sect. 2) we find that five sources are required. We re- 



xpix 



Fig. 5. X-ray contours in the central area of NGC 6397 as 
observed with the ROSAT HRI in 1995, with the positions 
of the sources obtained in the best fits. • positions found 
with the best fit for four components in C in which all 
positions are left free (Model I), □ positions found with the 
best fit for five components in C, in which three sources 
have distances fixed with respect to one other source, at 
the distances of ID 1, 3, and 6 to ID 2 (Model IV). Other 
symbols as in Fig.^. Modell is marginally better at the 
2-sigma level. 



do not detect source A of Cool et al. (1993) in the 1995 
observation. 

Cool et al. (1995b) resolved source X4/C into three 
components Cl-3. From a list of Ha emission and/or ul- 
traviolet excess objects (their Table 2), they suggest iden- 
tifications of IDl with C2, ID 2 for C3 and ID 3 for CI. 
Comparing the positions of the sources in Model I we find 
that the positions of X 4b and X 4c are compatible with 
those of ID 3 and ID 1, respectively; we thus identify X4b 
with CI and X4c with C2. X4d is a new source. (The 
offset required to match these positions from Table ^ with 
those given by Cool et al. (1995b) is slightly larger than 
our claimed accuracy of ^ 1"; the remaining difference 
may be explained by an offset between the Guide Star 
Catalogue coordinate system and the more accurate Hip- 
parcos coordinate system.) The position of X4a is not 
compatible with that of ID 2. The reason for this may 
be seen in Fig.^: the two brightest components of source 
X4/C have a smaller difference in right ascension than 
ID 2 and ID 3. If ID 2 is the correct identification for C3, 
we conclude that X4a is not identical to C3. 
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Table 4. Results of fitting four models to the three data 
sets of NGC 6397. The Table lists the number n of fit- 
ted parameters, and the difference AL with respect to the 
best model for a given data set. In all fits, the fluxes of 
all sources are fitted parameters. For the 1995 observation 
Model I has four sources with free positions. Model IV has 
six sources of which three have free positions and three 
have fixed positions relative to ID 3, corresponding to the 
offsets of ID 1, ID 2, and ID 6 with respect to ID 3. Mod- 
els II and III are as Model IV, after removing ID 6 and 
ID 2, respectively. For each Model, the same positions as 
in the best fit for 1995 are used for the 1991 and 1992 
data. 
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IV (B-f 1-flD 1236) 


13 
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1.9 



sources must match the distances between the proposed 
optical counterparts, which are accurately known from the 
HST observations. In Model II we fit five sources to the 
X-ray data of the center of NGC 6397, of which three are 
forced to be at fixed relative positions, corresponding to 
the distances between IDl, ID 2 and ID 3. Model II thus 
has four fitted parameters less than Model I. The InL of 
Model II is 26 higher than that of Model I, i.e. it is a sig- 
nificantly (4-sigma) worse fit. This confirms that X4a is 
not ID 2. In Model III we assume that ID 6 of Table 2 in 
Cool et al. (1995b) rather than ID 2 is the counterpart of 
C3, and fix the distances between the sources accordingly. 
This fit has the same In L as Model II, and thus also is sig- 
nificantly worse than the fit of Model I. Again, the reason 
for the bad fit is the mismatch in the difference in right 
ascensions of the two brightest X-ray sources with that 
between the proposed counterparts: X4a is not ID 6. 

We note that the best position of X 4a is between ID 2 
and ID 6, and in Model IV we fit six sources, of which four 
are forced to be at the relative distances of ID 1-3 and 
ID 6. Model IV thus has three fitted parameters less than 
Model I. Its In L is 6 higher than that of Model I, i.e. it is 
marginally worse at less than 2-sigma. The parameters of 
the six sources of this model are also given in Table |. It 
is seen that the positions of X4b, X4c and X4d are the 
same (within the error) in Model IV as in Model I. 

Thus, we have two acceptable models. In both mod- 
els we confirm the possible identifications of ID 3 with CI 
(=X4b) and of IDl with C2 (=X4c), and we find one 
new source (X4d). In Model I the remaining flux is as- 
cribed to one source (X4a) which is not identical to ID 2. 
In Model IV the remaining flux is ascribed to two sources, 
one of which is ID2/C3 and one is a second new source, 



Table 5. Countrates (counts ksec~^) assigned to the cen- 
tral sources in Models I and IV in the fits to the observa- 
tions of 1991 and 1992. Numbers in parentheses indicate 
the errors in the last digit. For 1995 see Table^. 
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4-3. The earlier observations 

The standard analysis detects X2, X5, X16, X6 and X8 
in both the 1991 and the 1992 data of NGC 6397, and X 19 
in the 1992 data, all at countrates compatible with those of 
1995. It also detects sources X 13/B and X4/C in the 1991 
data and in the 1992 data, labelling both as extended. The 
number of photons in sources B and C is rather small in 
these short observations. To limit the number of parame- 
ters in the fits to the central sources we demand that the 
distance between the fitted central sources in each model 
is the same as in the best fit to the 1995 data, but allow the 
fluxes to be different. The corresponding reductions in the 
number of fitted parameters for each model are indicated 
in Table|. 

We thus fit four models to each data set. For each year, 
the best model is set at AL = 0, and the quality of the 
other models for that year is determined with respect to 
this model. The results of our fitting are shown in Table ^ 
For the 1991 data, the models with five sources are com- 
parable in quality, and the six-source model is not signif- 
icantly better. For the 1992 data. Model III is marginally 
better (2 sigma) than Model I and significantly (3 sigma) 
better than Model II, whereas Model IV is of similar qual- 
ity. 

The fits to the earlier data confirm the conclusions that 
we draw on the basis of the observation of the long obser- 
vation of 1995. Modell in which source X4/C is separated 
into four components at free positions is acceptable for all 
three observations. Model II in which X4/C is separated 
into four components at fixed relative distances of ID 1- 
3, is not acceptable for the 1992 data. Model IV in which 
X4/C is separated into five components, four of which cor- 
respond to ID 1, ID 2, ID 3 and ID 6, also is acceptable for 
all observations. ID 2 is not required in 1992, and ID 6 is 
not required in 1991. The latter fact explains why ID 6 is 
not present in the analysis by Cool et al. (1993) of the 1991 
data. These conclusions are confirmed by the countrates 
that Model IV ascribes to the different sources, listed in 
Table I 



F. Verbunt & H.M. Johnston: Dim X-ray sources in four globular clusters 



9 



1992 observation (after shifting the 1992 observation by 
3'.'5, 1'.'O in Q, S; compatible with the shift as determined 
by Cool et al.). We fit Model I to the added imago, and 
compare it with the fit in which a source is added to 
model A. Wo find AL = 10, which implies that source A is 
marginally significant at ~ 2.5a. The position (+2!4, —14" 
with respect to source B) and countrate (0.4 counts/ks) 
that we find for source A are compatible with those given 
by Cool et al. (1993). 

4-4- Sources not related to the cluster 

HIP 86569 is a Kl IV/V star with V = 9.44, B-V = 0.90, 
and a parallax of 0.0167(18)". Hipparcos discovered that 
this star is a close binary (separation 0.19") of stars with 
Hipparcos magnitudes H = 10.17(9) and 10.47(11), re- 
spectively. At a distance of 60 pc the observed ROSAT 
HRI countrate converts to an X-ray luminosity in the 
0.5-2.5 keV band of Lx — 4 x lO^^ergs"-^ (for assumed 
1.4 keV bremsstrahlung with no absorption). This is sim- 
ilar to the X-ray luminosities of single KV stars detected 
in the ROSAT All Sky Survey, such as HD 17925 (KIV) 
which has ~ 1.2 x lO^^ ergs"^ (Hiinsch et al. 1998). 

4-5. Discussion 

The core of NGC 6397 contains at least four X-ray sources 
detected with ROSAT, and possibly five. 1 ctskscc"^ for a 
source at the distance and with the absorption column of 
NGC 6397, for an assumed 0.6 keV bremsstrahlung spec- 
trum corresponds to a luminosity in the 0.5-2.5 keV band 
of 2.2 X 10'^^ ergs~^. The faintest source we detect, X4c, 
is at this level. The brightest source is X 13/B, at a lumi- 
nosity of about 6.8 x lO'^^crgs^^. These luminosities are 
at the bright end of the luminosity distribution for cat- 
aclysmic variables, such as the large sample investigated 
with ROSAT (Verbunt et al. 1997), as expected for an 
X-ray selected sample. 

Of these sources, X13 and X4b have the same flux 
level in all thrc;e observations. Source X4c is fainter in 
1995. The identifications of X4b and X4c with cata- 
clysmic variables ID 3 and ID 1 remains probable, as does 
the argument by Edmonds et al. (1999) that these systems 
are DQ Her type systems. The distance between X4b and 
X4c in Model I is marginally less than the distance be- 
tween the proposed counterparts; it is tempting to specu- 
late that this is due to a small X-ray flux contribution of 
a fourth cataclysmic variable ('CV4') identified by Cool 
et al. (1998) and confirmed by Edmonds et al. (1999). 

The new source X4d has been detected in the 1995 
observation because of the longer exposure; it may, but 
need not, be brighter in the 1995 observation than in 1991 
and 1992. 

If the remaining flux is assigned to one source X4a, 



two sources X4a and X4e, the flux of X4a may still be 
constant, and X4a may be identified with the probably 
DQ Her type cataclysmic variable ID 2. In this case, the 
fliix of X4e has increased between 1991 and 1995. ID 6 
was reported to vary by 1.1 magnitude in five hours by 
De Marchi & Paresce (1994), but was constant in a ten 
hour observation by Cool et al. (1998). It is suggested by 
Edmonds et al. (1999) that ID 6 is a undermassive helium 
white dwarf, probably in a binary. If it is a single helium 
white dwarf, it cannot be a variable X-ray source; if it is 
in a binary with a recycled radio pulsar, it also is unlikely 
to be a variable X-ray source; if it is in a binary with an- 
other white dwarf, then optical and X-ray variability can 
be due to variable mass transfer from that other white 
dwarf. However, it is also possible that not ID 6, but a 
nearby hitherto unidentified star in NGC 6397 is the X- 
ray source X4e. 

Whether source X4a alone, or source X4a and X4e 
are present in the core of NGC 6397, and in the latter case 
whether X 4e is identical to ID 6 requires a better spatial 
resolution for the X-ray observations than provided by 
ROSAT. 

5. NGC 6752 

Two sources have been detected in the core and two near 
the core of NGC 6752 in a ROSAT HRI observation ob- 
tained in 1992 (Grindlay 1993); close to one of the core 
sources, two candidate cataclysmic variables have been 
identified on the basis of H a emission and variability on 
(presumably orbital) periods of 5.1 and 3.7 hrs (Bailyn et 
al. 1996). 

5.1. Data analysis and source list 

Two more observations of NGC 6752 have been obtained 
by us. Because the three observations have comparable 
length, we add them into a combined image which we 
analyse and use as a reference for the individual obser- 
vations. To add the three observations we use the method 
outlined by Verbunt & Hasinger (1998), as follows. First 
we correct the data for each observation separately for 
the changed pixel size (see Sect. 2), analyse the result- 
ing images and determine the offsets between sources de- 
tected in separate observations. Averaging these offsets 
we find (on the basis of sources X3, X13, X4, and X6) 
that the X-ray coordinates of the 1996 observation have 
to be shifted by da = — 0'.'7 ± 0'.'7 in right ascension 
and dS = +3'.'6 ± 0'.'7 in declination to be brought in 
line with the 1992 data. Similarly, the 1995 data (on the 
basis of the same sources plus X14) must be shifted by 
da ^ +2'.'4: ± 0"7, dS = +1'.'2 ± 0'.'7. We apply these cor- 
rections to the pixel coordinates of the photons, and then 
add the three images into a combined image, which is anal- 
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Table 6. X-ray sources detected in NGC 6752 {Ay = 0.12, 
d = 4.2 kpc, Djorgovski 1993) with the ROSAT HRI, for 
the standard analysis of the whole field, and separately 
for the multi-source analysis of the central area. Numbers 
up to 15 are sources from Johnston ct al. (1994), higher 
numbers are new. All X-ray positions have been corrected 
for boresight. Identifications with letter on the right refer 
to Grindlay (1993). The positions of the center of the clus- 
ter (GC, Djorgovski & Meylan 1993), it core radius and 
half-mass radius (Trager et al. 1993) and the positions of 
some optical objects discussed in the text are also listed. 
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We identify two sources with stars with accurate po- 
sitions: X19 is close to TYC 90712281 (CD-60 7128), a 
star with V = 9.99 and B = 10.48, and Xll is close to 
HIP 94235/HD 178085, a GOV star with V ^ 8.38 and 
B = 9.00. The latter identification was suggested already 
by Johnston et al. (1994) on the basis of the PSPC obser- 
vation. The chance probability of finding a counterpart at 
these optical brightnesses is small, and we consider both 
identifications secure. The X-ray position of Xll is af- 
fected by its proximity to the edge of the HRI detector. 
For this reason we use X 19 to tie the X-ray to the optical 
coordinates. X 19 is not found in any of the three individ- 
ual observations, showing up only in the combined frame. 
It is a relatively weak source and its position accordingly 
has an error of I'.'S both in right ascension and in declina- 



100 -100 




xpix 



Fig. 6. X-ray contours in the central area of NGC 6752 as 
observed with the ROSAT HRI in the combined image of 
the 1992-1996 data. The image was smoothed with a 2-d 
(T~3" Gaussian. The detected sources are indicated with 
their numbers in Table||. The circle gives the half-mass 
radius of the cluster. The lower and left axes give pixel 
numbers for the ROSAT HRI detector, the upper and right 
axes right ascension and declination with respect to the 
cluster center. 

to all positions of the X-ray sources; the resulting positions 
are listed in Table |[ (The remaining offset between Xll 
and HIP 94235 is within the nominal error for the right as- 
cension, and within 2-sigma for the declination: note that 
the error is composed of the statistical uncertainties in the 
positions of both X 19 and X 11.) 

The detection limit in the total observation is about 
0.7 X 10^^'* ergcm~^s~^. An area with radius 12.'5 in the 
ROSAT Deep Survey contains 25 sources brighter than 
this limit; we thus expect to find 0.6 in the region within 
the half-mass radius of NGC 6752, rh — 2' . The sources in 
the core thus probably belong to the cluster, and possibly 
the two sources X6/A and X14/D as well. 

5.2. Sources in the center of the cluster 

Analysing the central source with the method described 
in Sect. 2, we find four significant sources (the increase in 
InL is 29 both for the third and for the fourth source). 
This adds two sources to the two already described by 
Grindlay (1993). The position and fluxes of these sources 
are listed in Table IbI Fig. ffl shows the positions and X- 
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Fig. 7. X-ray contours in the central area of NGC 6752 as 
observed with the ROSAT HRI in the combined data of 
1992-1996, with the positions of the sources obtained with 
the multi-source fit (□). The positions of two candidate 
cataclysmic variables found by Bailyn et al. (1996) are 
indicated with 

by Bailyn et al. (1996). The southern cataclysmic variable 
('star r) is at 3.8 ±2.3" from X7a, and therefore remains 
a possible counterpart (assuming an error of 1" for the 
optical position, and taking into account the 2" error of 
X19). 

We have analysed the separate observations, keeping 
the position of the four central sources fixed at those of 
the co-added image (as listed in Table ^, but allowing 
their fluxes to vary. We do not find significant evidence for 
variation; due to the limited statistics we cannot exclude 
variations by a factor two. 

A countrate of 1 ctsksec"^ corresponds to a lumi- 
nosity between 0.5 and 2.5 keV of 4.6 x lO'^^ergs"^ at 
the distance of NGC 6752 and for an assumed 0.6 keV 
bremsstrahlung spectrum. Thus, X7a and X7b have X- 
ray luminosities of about 7.5 x lO^^ergs"^, and X21 
and X 22 about a quarter of this. X 6 and X 14, the 
two sources outside the cluster core have luminosities of 
4.6 X lO'^^ergs"^ and 3.2 x lO'^-^ ergs"-^, respectively, if 
they are in NGC 6752. 

5.3. Sources not related to the cluster 

The spectral type of TYC 9071 228 1 is not known; on the 
basis of its magnitude and colour (^ = 9.99,-6 — ^ = 



1.4 keV bremsstrahlung spectrum, the countrate of X19 
converts to an X-ray luminosity in the 0.5-2.5 keV band of 
~ 3 X lO^^ergs"^, a reasonable value for a late F main- 
sequence star (see e.g. the fist of ROSAT detections of 
bright stars by Hiinsch et al. 1998). The Tycho Catalogue 
marks this star with 'unresolved duplicity', with visual 
magnitude varying between 9.51 and 10.88. 

HIP 94235 has a significant parallax which puts it at 
57 pc. Its countrate converts to an X-ray luminosity at 
that distance of about 2 x 10^^ erg s^^, a normal X-ray 
luminosity for a GOV star. 

Comparison of the ROSAT image with the USN0-A2 
Catalogue gives a candidate identification for X 16, at a 
distance of 2'.'2, see Table^. No other sources outside the 
cluster have been identified by us. 

We have analysed the three separate HRI observations, 
and find no evidence for variablility, except for X 15, which 
in March 1992 had an X-ray fiux about half of that ob- 
served in March 1995 and April 1996. 

6. Liller 1 

The globular cluster Liller 1 is a highly reddened cluster 
near the galactic center [Ay ~ 9.5, d — 8.6 kpc, Frogel et 
al. 1995). It probably has undergone core collapse (Djor- 
govski 1993). Liller 1 harbours the Rapid Burster, a highly 
unusual recurrent transient. When discovered in 1977 the 
source emitted short (^ 5 s) bursts of X-rays every ^ 10 s; 
in some later observations, e.g. Aug 1985, it emitted bursts 
of ~ 500 s separated by 1500-4000 s; and it has also been 
observed as a steady source. The bursts are interpreted 
as accretion events. In addition to these, thermonuclear 
bursts have also been detected, identifying the accreting 
star as a neutron star. A review of this remarkable source 
is given by Lewin et al. (1995). A low-luminosity X-ray 
source near Liller 1 is tentatively identified as the quies- 
cent (low-state) counterpart of the Rapid Burster (Asai et 
al. 1996). 

No source is detected in the cluster in our ROSAT 
HRI observation of the globular cluster Liller 1. Near the 
cluster center, no circle with radius of 5" contains more 
than 4 photons. For an expected number of 10 photons, 
the probability of getting 4 or fewer photons is less than 
4%. We thus take 10 as the 2-cr upper limit to the num- 
ber of photons, which with the effective exposure time is 
converted to an upper limit of 0.6 ctsksec"^. 

Asai et al. (1996) report the detection on 1993 Aug 27 
with ASCA of a source near Liller 1. For a powerlaw with 
photon index 2, absorbed by a column A^h — lO^^cm"^, 
this source has an unabsorbed fiux in the 2-10 keV band 
of 2.5lQg X 10~^^ ergcm~^s~^. For this spectrum our up- 
per limit in the ROSAT HRI corresponds to a flux of 
1.4 X 10"-^'^ ergcm^^s^^, slightly lower than the ASCA 
detection. 
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Table 7. Positions of the center of Liller 1 (GC, Picard & 
Johnston 1995) its core radius (Trager et al. 1993) and the 
positions of X-ray sources detected near it, viz. the Rapid 
Burster (RB, Hertz & Grindlay 1983), a dim source de- 
tected with ASCA (A, Asai et al. 1996) and a dim source 
detected with the ROSAT HRI (R, this paper). The po- 
sition of the O star HD 317889 is also given. The final 
columns gives the errors in the positions. 





a (2000) 


5 (2000) 


San 


55(") 


GC 


17 33 24.5 


-33 23 20.4 


re = 4" 




RB 


17 33 24.0 


-33 23 16.2 


2 


2 


A 


17 33 19.0 


-33 23 


90 


30 


A" 


17 33 24.2 


-33 23 6.5 


20 


20 


R 


17 33 4.76 


-33 23 27.2 


5 


5 


HD 


17 33 5.02 


-33 23 28.4 







"new determination by Asai, 1999, private communication 

o) Cen o • . . o 

NGC6397 o .... 
NGC6752 . . o o .. 
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l°g Lo.5_2.5 kev (erg/s) 

Fig. 8. X-ray luminosities of dim sources in four globu- 
lar clusters. Sources in and outside the core are shown as 
• and o, respectively. Data points are from Table |[ and 
for 47 Tuc from Verbunt & Hasinger (1998) slightly mod- 
ified for an assumed 0.6 keV bremsstrahlung spectrum. In 
all cases the detection limits in and outside the core are 
close to the lowest detected luminosities in and outside 
the cores. 

error in the position of this source is about 1"; the actual 
error is dominated by the error in the bore sight correction, 
which is about 5". The ROSAT source is not compatible 
with the center of Liller 1, and also not compatible with 
the position of the Rapid Burster as determined with Ein- 
stein (see Table|). The position of the ROSAT source co- 
incides within the bore sight uncertainty with the 04 III(f ) 
star HD 317889 (Vijapurkar & Drilling 1993). The star is 
in the Tycho Catalogue as TYC 7380976 1. From the ob- 
served magnitude and colours {V = 10.12, B — V ~ 0.92, 
U—B = —0.23, Drilling 1991) we estimate a reddening and 
distance of E{B - F) ~ 1.2 and d ~ 3kpc for HD 317889. 
The observed ROSAT HRI countrate is as expected for 
such a star, according to the general correlation between 
bolometric luminosity and X-ray luminosity of O stars: 
~ IQ-'^Lb (e.g. Kudritzki et al. 1996). (HD 317888 is 
within 1" of the 04 star; we have not been able to find 
more information on this star.) 



did detect the Rapid Burster in quiescence, or another low- 
luminosity source in Liller 1; and that ROSAT observed 
when this source had a lower flux level. In fact, varia- 
tion of transients in their quiescent state is common (e.g. 
Campana et al. 1997). The star detected with ROSAT 
in this case is not detected with ASCA, presumably be- 
cause its spectrum is too soft. The second interpretation 
is that ASCA in fact detected the star also detected with 
ROSAT, and not the quiescent counterpart of the Rapid 
Burster. The position of the ROSAT source is marginally 
compatible with that of the ASCA source; its countrate is 
exactly that predicted on the basis of the ASCA source. 
Dr. Asai has kindly communicated a new determination of 
the position of the X-ray source detected by ASCA, using 
new calibrations to improve the accuracy. This position, 
listed in Table |^, excludes the ROSAT source as a pos- 
sible counterpart, and thus confirms that ASCA indeed 
detected a source in the cluster. 

7. Summary and discussion 

In the three low-reddened clusters uj Cen, NGC 6397 and 
NGC6752 we have detected a total of 17 dim X-ray 
sources, of which 5 are well outside the core. The X-ray 
luminosities of these sources are listed in Table^, and plot- 
ted in Fig.^. The interpretation of Fig.^ must be made 
with some care. First, sources outside the core may not 
belong to the cluster; the faintest core source in uj Cen 
may be a fore- or background source. Second, the con- 
version of observed countrate to luminosity depends on 
the assumed spectrum, and from PSPC observations we 
know that different sources have different spectral param- 
eters (Johnston et al. 1994). For example, the 0.6 keV 
black body spectrum used for the sources in lo Cen gives 
a 40% higher flux for the same countrate than an as- 
sumed 0.6 keV bremsstrahlung spectrum would give. The 
bremsstrahlung spectrum is used for the three other clus- 
ters. Third, the detection limits in NGC 6397, NGC 6752 
and 47 Tuc are higher in the cores, where the point spread 
functions of sources overlap, than outside the core. Such 
a difference is not present in uj Cen. Fourth, we show the 
average luminosity, and several sources are known to be 
variable. 

With these points in mind, we note from Fig.^ that 
in all clusters except possibly lo Cen the most luminous 
sources appear to be in the cluster core. The main dif- 
ference between uj Cen and the other clusters is that the 
collision frequency in oj Cen is so low that one expects no 
low-mass X-ray binaries in it, and that most cataclysmic 
variables in it will be evolved from primordial binaries 
(Verbunt & Meylan 1988, Davies 1997). In addition, the 
mass segregation in this cluster is very low. Thus in oj Cen 
there is no marked difference between the core and the 
regions outside the core. 
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Table 8. X-ray luminosities in ergs^^ in the 0.5-2.5 keV 
band of the dim X-ray sources in globular clusters de- 
scribed in this paper. For sources in lo Cen we assume 
a 0.6 keV black body spectrum; for those in NGC 6397 
and NGC 6752 a 0.6 keV bremsstrahlung spectrum. For 
the same countrate, the blackbody spectrum corresponds 
to a flux higher by about 40% than the bremsstrahlung 
spectrum. 



uj Cen 


NGC 6397 


NGC 6752 


X logLx 


X 


log Lx 


X 


log Lx 


core 




core 




core 


9a 32.14 


13 


31.84 


7a 


31.86 


9b 32.22 


4a 


31.73 


7b 


31.89 


20 31.96 


4b 


31.75 


21 


31.35 


outside 


4c 


31.43 


22 


31.26 


7 32.33 


4d 


31.53 




outside 


21 31.88 




outside 


6 


31.66 




12 


30.95 


14 


31.50 



will detect more sources. In the cores of NGC 6397 and 
NGC 6752 the detection of more source will also require 
better imaging, so that the faint sources can be detected 
against the brighter ones. We do not detect a difference be- 
tween the luminosities of sources in the collapsed globular 
cluster NGC 6397 and the much less concentrated globular 
cluster NGC 6752. On the other hand, the highly concen- 
trated cluster 47 Tuc contains three sources which are an 
order of magnitude brighter than the brightest sources in 
NGC 6397 and NGC 6752. 

Viable optical counterparts have been suggested for 
only five among the 26 sources shown in Fig. ^, all of them 
probable cataclysmic variables. We compare the ratio of 
X-ray flux to optical flux of these sources with the ratios 
measured for cataclysmic variables and for RS CVn sys- 
tems in the Galactic Disk in Fig.||. It is seen that the sug- 
gested optical counterparts for the sources in NGC 6397 
and NGC 6752 lead to ratios which are compatible with 
those of cataclysmic variables, whereas those in 47 Tuc 
are too bright in X-rays, in agreement with Fig.^. If these 
sources are indeed cataclysmic variables, their excessive 
X-ray luminosity needs to be explained; alternatively, the 
suggested identifications may be chance coincidences (as 
discussed by Verbunt & Hasinger 1998). All suggested 
counterparts lead to higher X-ray to optical flux ratios 
than those of RS CVn binaries. 

The accurate positions that we determine for individ- 
ual sources are valid for separately detected sources in 
particular. In the case of overlapping sources, we do not 
have unique solutions. Thus, in the core of NGC 6397 fits 
with 5 and 6 sources are both acceptable, at similar qual- 
ity; and we cannot exclude that more sources contribute 
to the observed fiux, which would invalidate our derived 
positions. 
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Fig. 9. X-ray countrates of the dim sources in globular 
clusters as a function of their visual magnitude, compared 
with the ROSAT PSPC countrates and visual magnitudes 
of various types of cataclysmic variables (data from Ver- 
bunt et al. 1997; filled symbols represent systems first dis- 
covered in X-rays and only subsequently identified as cat- 
aclysmic variables, i.e. X-ray selected systems) and with 
RS CVn systems (data from Dempsey et al. 1993) respec- 
tively. PSPC countrates of the dim cluster sources have 
been computed for an assumed 0.6 keV bremsstrahlung 
spectrum, corrected for absorption, from the observed HRI 
countrates. Visual magnitudes are also corrected for ab- 
sorption. T indicate sources in 47 Tuc (X9 and X19, V as 
estimated by Verbunt & Hasinger 1998), A in NGC 6397 
(X4b and X4c, V from Cool et al. 1998), B in NGC 6752 
(X7a, V from Bailyn et al. 1996). The dotted lines indi- 
cate a constant ratio of X-ray to optical flux. 

cause a three-body interaction (i.e. a close encounter of a 
binary with a single star) in the core has expelled the bi- 
nary from the core (e.g. Hut et al. 1992). In the latter case 
the binary is expected to be eccentric immediately after 
being expelled; tidal forces may in time circularize the or- 
bit again. Such binaries are only a minority of the overall 
binary population of a cluster; however. X-ray observa- 
tions may preferably select such binaries if tidal forces act 
in them. Since sources away from the core can be fore- 
or background sources, optical identification of them is 
required to settle whether they belong to the cluster or 
not. Our accurate positions should help in finding such 
counterparts. 

Acknowledgements. We have made use of the ROSAT Data 
Archive of the Max Planck Institut fiir extraterrestrische 
Physik at Garching; of the SIMBAD database operated at 
Centre de Donnees astronomiques in Strasbourg; and of the 
Digitized Sky Surveys. For Hipparcos and Tycho data the CD- 
ROM Celestia. orovided bv ESA. was verv helDful. We further 



14 



F. Verbunt & H.M. Johnston: Dim X-ray sources in four globular clusters 



the Digitized Sky Surveys, and Dr. Asai for communication 
of the re-determination of the position of the ASCA source in 
Lillerl. 

References 

Asai, K., Dotani, T., Kunieda, H., Kawai, N. 1996, Publ. Astr. 
Soc. Japan, 48, L27 

Bailyn, C, Rubenstein, E., Slavin, S., et al. 1996, ApJ (Let- 
ters), 473, L31 

Berghofer, T., Schmitt, J., Hiinsch, M. 1999, A&A, 342, L17 
Campana, S., Mereghetti, S., Stella, L., Colpi, M. 1997, A&A, 

324, 941 
Cash, W. 1979, ApJ, 228, 939 

Cool, A., Grindlay, J., Bailyn, C, Callanan, P., Hertz, P. 1995a, 

ApJ, 438, 719 

Cool, A., Grindlay J., Cohn, H., Lugger, P., Bailyn, C. 1998, 

ApJ (Letters), 508, L75 
Cool, A., Grindlay, J., Cohn, H., Lugger, P., Slavin, S. 1995b, 

ApJ, 439, 695 

Cool, A., Grindlay, J., Krockenberger, M., Bailyn, C. 1993, 

ApJ (Letters), 410, L103 
Cruddace, R., Hasinger, G., Schmitt, J. 1988, in F. Murtagh, 

A. Heck (eds.). Astronomy from large databases, p. 177 
Danner, R., Kulkarni, S., Thorsett, S. 1994, ApJ (Letters), 436, 

L153 

David, L., Harnden, Jr, F., Kearns, K., Zombeck, M. 1995, The 
ROSAT High Resolution Imager (HRI), Technical report, 
U.S. ROSAT Science Data Center/SAO 
Davies, M. 1997, MNRAS, 288, 117 
de Marchi, G., Paresce, F. 1994, A&A, 281, L13 
Dempsey, R., Linsky, J., Fleming, T., Schmitt, J. 1993, ApJS, 
86, 599 

Djorgovski, S. 1993, in S. Djorgovski, G. Meylan (eds.), Struc- 
ture and Dynamics of Globular Clusters, Vol. 50 of ASP 
Conference Series, ASP, p. 373 

Djorgovski, S., Meylan, G. 1993, in S. Djorgovski, G. Mey- 
lan (eds.). Structure and Dynamics of Globular Clusters, 
Vol. 50 of ASP Conference Series, ASP, p. 325 

Drilling, J. 1991, ApJS, 76, 1033 

Edmonds, P., Grindlay, J., Cool, A., et al. 1999, ApJ, 516, 250 
ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP 
1200 

Frogel, J., Kuchinski, L., Tiede, G. 1995, AJ, 109, 1154 
Grindlay, J. 1993, in S. Djorgovski, G. Meylan (eds.). Struc- 
ture and Dynamics of Globular Clusters, Vol. 50 of ASP 
Conference Series, ASP, p. 285 
Grindlay, J., Cool, A., Callanan, P., et al. 1995, ApJ, 455, L47 
Hakala, P., Charles, P., Johnston, H., Verbunt, F. 1997, MN- 
RAS, 285, 693 

Hasinger, G., Burg, R., Giacconi, R., et al. 1998, A&A, 329, 
482 

Hertz, P., Grindlay, J. 1983, ApJ, 275, 105 

H0g, E., Biissgen, G., Bastian, U., et al. 1997, A&A, 323, L57 

H0g, E., Kuzmin, A., Bastian, U., et al. 1998, A&A, 335, L65 

Hiinsch, M., Schmitt, J., Voges, W. 1998, A&AS, 132, 155 

Johnston, H., Verbunt, F. 1996, A&A, 312, 80 

Johnston, H., Verbunt, F., Hasineer, G. 1994, A&A, 289, 763 



Kaluzny, J., Kubiak, M., Szymanski, M., et al. 1997, A&AS, 
122, 471 

Koch-Miramond, L., Auriere, M. 1987, A&A, 183, 1 

Kudritzki, R. P., Palsa, R., Feldmcier. A.. Puis, J., Pauldrach, 
A. W. A. 1996, in H. Zimmermann, J. Triimper, H. Yorke 
(eds.), Roentgenstrahlung from the Universe, MPE Report 
263, p. 9 

Lewin, W., van Paradijs, J., Taam, R. 1995, in W. Lewin, 

J. van Paradijs, E. van den Heuvel (eds.), X-ray binaries, 

Cambridge U.P., Cambridge, p. 175 
Margon, B., Bolte, M. 1987, ApJ, 321, L61 
Mattox, J., Bertsch, D., Chiang, J., et al. 1996, ApJ, 461, 396 
McGale, P., Pye, J., Hodgkin, S. 1996, MNRAS, 280, 627 
Monet, D., Bird, A., Canzian, B., et al. 1998, The PMM 

USNO-A2.0 Catalog, US Naval Observatory 
Ferryman, M., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 

323, L49 

Picard, A., Johnston, H. 1995, A&AS, 112, 89 

Trager, S., Djorgovski, S., King, I. 1993, in S. Djorgovski, G. 
Meylan (eds.). Structure and Dynamics of Globular Clus- 
ters, Vol. 50 of ASP Conference Series, ASP, p. 347 

Triimper, J., Hasinger, G., Aschenbach, B., et al. 1991, Nat, 
349, 579 

Verbunt, F., Bunk, W., Hasinger, G., Johnston, H. 1995, A&A, 
300, 732 

Verbunt, F., Bunk, W., Ritter, H., PfefFermann, E. 1997, A&A, 

327, 602 

Verbunt, F., Hasinger, G. 1998, A&A, 336, 895 

Verbunt, F., Meylan, G. 1988, A&A, 203, 297 

Verbunt, F., Shafer, R., Jansen, F., Arnaud, K., van Paradijs, 
J. 1988, A&A, 168, 169 

Vijapurkar, J., Drilling, J. 1993, ApJS, 89, 293 

Zimmermann, H., Becker, W., Belloni, T., Dobereiner, S., Izzo, 
C, Kahabka, P., Schwentker, O. 1996, EXSAS User's 
Guide: Extended scientific analysis system to evaluate data 
from the astronomical X-ray satellite ROSAT, Edition 5, 
Technical report, MPE 



